Universal features of fluctuations 
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Universal scaling laws of fluctuations (the A-scaling laws) can be derived for equilibrium and off- 
equilibrium systems when combined with the finite-size scaling analysis. In any system in which 
the second-order critical behavior can be identified, the relation between order parameter, criticality 
and scaling law of fiuctuations has been established and the relation between the scaling function 
and the critical exponents has been found. 



I. INTRODUCTION 

Recent years have witnessed an intense experimental 
and theoretical activity in search of scale-invariance and 
fractality in multihadron production processes |l||. The 
creation of soft hadrons in these processes is related to the 
strong-coupling, long-distance regime of Quantum Chro- 
modynamics (QCD) whose exploration remains the sa- 
cred grail of the high-energy particle physics. Through a 
multitude of experimental studies of particle density fluc- 
tuations and particle correlations, new informations have 
been gathered which prompted the suggestion to investi- 
gate the pattern of multiplicity fluctuations in the small 
domains of phase-space looking for the scale-invariance 
and fractality pj. New measurements has triggered re- 
vival of interest in the theory of fluctuations in complex 
systems. 

In this work we shall discuss universal features of fluc- 
tuations of physical quantities in a Z?-dimensional, A'^- 
body system, where N is essentially finite. We are par- 
ticularly interested in the self-similar systems, such as 
the fractal objects or the thermodynamic systems at the 
second-order phase transition. By the self-similarity we 
mean that the characteristic length ~ (7V*)^/^ {N* is 
the characteristic size), which could be associated with 
the disappearance of fluctuations cannot be defined. The 
progress in this domain is associated with the develop- 
ment of the renormalization group methods. Actually, 
there is a close connection of the renormalization group 
ideas and the limit theorems in the probability theory. 
This basic idea stems from the works of Bleher, Sinai ||] 
and Jona-Lasinio Q , and consists of splitting the whole 
system into correlated blocks. Then the probability dis- 
tribution of the investigated physical quantity is calcu- 
lated for the blocks, and the renormalization procedure 
is expressing the distribution for large blocks in terms of 
the distributions for the smaller ones. These probability 
distributions dealing with larger and larger number of 
subunits are assumed to be given by probabilistic limit 



laws and the parameters of these distributions, properly 
renormalized at each step, are expected to tend to some 
finite limit. When this scheme is realized, the renor- 
malization procedures gives the universal behavior of the 
infinite system. 

Following this approach, we shall discuss how the uni- 
versal scaling laws of fluctuations of different observables 
appear in finite systems. In particular, we shall consider 
the order parameter fluctuations in any equilibrium or 
non-equilibrium system in which the second-order criti- 
cal behavior can be identified. These considerations will 
provide an understanding of the relation between the or- 
der parameter, the criticality and the scaling law of fluc- 
tuations. The notion of the relevant observable will also 
appear from this discussion. 



II. BROAD DISTRIBUTIONS AND 
LONG-RANGE CORRELATIONS 

A. A basic example of the probabilistic limit laws 

Most physically interesting 'macroscopic variables' are 
defined as sum of many 'local random variables'. To un- 
derstand such a 'macroscopic behavior' is to seek for some 
universal limit behavior. The simplest, non-trivial exam- 
ple is when the local variables are uncorrelated. Suppose 
that we make N trials and a certain event has finite prob- 
ability Po to occur at each trial. In this chain of trials, we 
are looking for the behavior of the probability distribu- 
tion PN[m\ to get exactly m events during N trials. In 
the following, m will be called the multiplicity. This kind 
of process could describe for example the successive emis- 
sion of particles : at each time interval r, one particle is 
emitted with a probability po from a parent body. The 
equation for the probability distribution of the number 
of emitted particles in the time interval nr : 



PN[m] = (1 -po)PN-i[m] +PoPN-i[m- 1] 



(1) 



tells that the multiplicity at stage N is equal to m if, 
either it was m at stage N — 1 and then the emission 
does not occur at the following trial, or it was ?7i — 1 at 
stage A^ — 1 and the emission occurs in the last trial. In 
particular, eq. (|]) shows that the average multiplicity is 
< m >— Npo- When the number of trials N becomes 
large, the probability distribution Pj^[m\ tends to a nor- 
mal law : 



lim < TO >-^/^ PN[m]- 

<m> — >oo 



Up^ cxp 



-6„ 



m~ < m > 



< m > 



1/2 



(2) 



with bp^ = 1/[2(1 — Po)]- This form of the normal 
law expresses a relation between the scaled distribution 
< TO >^/^ P/v[?7i], and the shifted scaled multiplicity pa- 
rameter (to— < m >)/ < TO >^". 

Let us now introduce the correlation between the two 
consecutive trials by prohibiting the successive emission. 
In this case, the probability distribution verifies : 



Fat [to] = (l-p)Pjv-iH +pPN^2[m~l] 



(3) 



It is a simple exercise to obtain the solution of eq. (g) : 
the leading order of the average multiplicity is < ttt. >^ 
np/{l+p) H and the limit distribution is given again by a 
normal law (||) with bp^ — 1/[2(1 — Po)(l — '2po)]- Hence, 
the short-range correlations in (ph reduce the variance 
of the probability distribution by a factor 1 — 2po. It 
is important to realize a fundamental difference between 
the limit law, which is universal in the sense that most 
details of the considered process are unessential, and the 
parameters of the limit laws, which strongly depend on 
the specific features of the process. 



B. Central limit and stable laws for non-correlated 
variables 

Let us consider the additive, 'global' quantity : Y/v = 
^j^^i Xn, expressed in terms of 'local' microscopic vari- 
ables Xn- The fundamental problem in the statistical 
description of composed objects is to find a limit dis- 
tribution of this variable, with N being the number of 
subunits of the system. This problem can be approached 
in two steps : (i) when the random variables Xn are all 
independent, then (ii) when they are correlated. The first 
step has been solved entirely |^. The second step, when 
the variables are correlated, is not yet completely solved 
but some precise results and tools are available (see Ref. 
M for recent references). 

In the following, we will work with the probability den- 
sity : f{x)dx — Prob[x < X < x + dx]^ its characteristic 
function : 



4>x{k) = 



'"/(x)da 



and the expansion of the characteristic function in a series 
of cumulants near fc = : 



ln0x(fc) ^ik < X > 



— (<X2 >- <X>2) + 



(4) 



The importance of the characteristic function comes from 
the convolution theorem which tells that, if Xi and X2 
are two independent random variables and 4>Xi (k) and 
4'X2 (k) their characteristic functions, then the character- 
istic function of the random variable Y = C1X1+C2X2 is : 
4>Y{k) = 4>Xi{cik)4>x2{c2k) (ci and C2 are real numbers). 
In other words, the density of the sum of two indepen- 
dent variables is the product of the convolution of the 
two probability densities. The key point to determine all 
possible limit distribution of the random additive vari- 
able is the 'stability problem' : What are all probability 
densities f{x), such that, if Xi and X2 are two indepen- 
dent random variables of a probability density f{x), and 
Ci and C2 are two real numbers, then Y — CiXi +C2X2 is 
a random variable with the same probability density f{x) 
? It turns out that the stable distributions are those 
for which the logarithm of their characteristic function is 
given by the formula : 



\TL(l){k) = i-ik ~ c|fc|^ ( 1 + il3—uj^{k) 



(5) 



with four real parameters : j,c,(3 and fi, where c > 0, 
\/3\ < 1 and < /i < 2. In addition : uj^{k) — tan(7r/i/2) 
if /i ^ 1, else uJi{k) = 21n(fc/7r). This four-parameter 
family of functions can be written in different forms. One 
of them is : 



P;.,/3(ci/'^x + 7) 



1 



"1/m 



X Re 



if /i 7^ 1, and 



exp 



(~ikx - (1 + i/3tan — )fc^] dk 



Pijj{cx + ^l) = X 

TTC 



2/3 



cxp(— fc) cos ( k{x ^ — — In k) ) dk 



if /i = 1. For |a;| -^ cx), P^^p ex |a;|-(i+^) (0 < // < 2). 
For < ^ < 1 and (3^0, there is an essential singularity 
when X -^ : P^_/3 cx exp[— a(/i,/3)a;^/*^"'^~''-']. For other 
values, Pfi^fi is finite at the origin. 

The parameters appearing in these formulas are not 
equally important. 7 corresponds to a shift of the whole 
distribution, while c is a scale factor. These two param- 
eters are related to a translation and a dilatation of the 
distribution and can always be put respectively equal to 
and 1 by a suitable choice of the coordinates. The param- 
eter 13 is related to the asymmetry of the distribution and 



P vanishes for an even distribution. The most important 
parameter is the characteristic exponent /i. It cannot 
be larger than 2 because P^,^ should be non-negative, 
and is larger than to ensure the convergence at the 
origin. The normal (gaussian) distribution corresponds 
to /i = 2, and this is the only case where the variance 
ax =< X^ > — < X >'^ is finite. In fact, using (^ and 
the form of the stable distributions, one may show that : 



</i< 1 

1 <M< 2 

fi = 2 



< X >= oo and ax — oo 

< X >< oo and ax — oo 

< X >< oo and ax < oo 



D. Stable distributions and the self-similarity 

The stochastic process Yx is called self-similar with 
the Hurst exponent H if both Y\x and X^Yn have the 
same distribution. If Y/v is given by (||) , (0) , respectively, 
and the probability distribution P[3^aa] converges to the 
stable law P^,/3, then the stochastic process Yx is self- 
similar with the Hurst exponent H = l//i, and < Yx >~ 
TVI/ai = ]\f9 ^ g jg fj^Q anomalous dimension related to the 
extensive variable Y/v- 



E. Concept of the A-scaling 



C. Asymptotically stable laws - the domains of 
partial attraction 

Let us recall that the stable distribution is invariant 
when added to itself. One would like to know what kind 
of distributions of the local random variable X„ are lead- 
ing to the stable distribution P^^p when N ^ oo. Let 
Xi,X2, ■■■,Xx denotes a sequence of independent num- 
bers distributed according to a density f{x), and let Y/v 
be a random variable : 



Yn=Xi 



X 



and 



Y, 



N 



Xi 



N 



X 



N 



if 0< fi<l 



N <X > 



(6) 



(7) 



if 1 < /^ < 2. Let us define : yx — Yx/Bx, with the 
parameter Bx chosen in such a way that the probability 
distribution of 3V converges to the stable distribution 
Pfj.,/3 when N tends to the infinity. When the distribution 
f{x) yields the limit density P^./j for 3^, then / is said to 
belong to the domain of attraction of P^,^. The ensemble 
of the probability densities f{x), such that : 



lim 



lim 



Cfix)dx + JZf{x')dx' 

Cf{x)dx + j:::f{x')dx' 



lZn^')dx' I -13 



Cfi^')d^' 



■f3 



constitutes the entire domain of attraction of Pfj,.f3. We 
see that only the tail of the distribution f{x) is signifi- 
cant. If f{x) ~ \x\-^^+i''> (0 < /x < 2) for |a;| -^ oo, 
then f{x) belongs to the domain of attraction of Pp,/3, 
i.e. the reduced global variable y^ has a limiting stable 
distribution of index /x when N —^ oo. Moreover, all pa- 
rameters Bx,"f,c,P are known explicitly. In particular, 
Bx -- iVi/^. 



The case : 1 < /i < 2, corresponds to an infinite vari- 
ance of microscopic variable X with finite mean value 
< X >. Let Mx = Xi + ... + Xx. In this case : 



< Mx >^< Xi > 



<Xx >=N <X > 



The preceding characterization of the attraction domains 
of the stable laws yields : 



Prob 



< 



Mx - N < X > 



B 



N 



<y + dy 



P^,f3iy)dy 



with Bx ^< Mx >^/^. Introducing the probability den- 
sity /mjv('7i) of the variable Adx, one obtains the limit 
distribution : 

<M.>V./,,„M.P^,,(!^_±|^) (8) 

, when the number N of subunits X„ tends to the infinity. 
Notice that P^^p appears now as the scaling Junction with 
the scaling index A = l//i, which can be extracted by 
plotting < Mat >i/^ /mjv("^) against a scaling variable 
(to— < Mx >)/ < Mx >^l^ . In this representation, 
different curves for systems or stochastic processes with 
different number of subunits X collapse into a single scal- 
ing curve which characterizes the stable law. The scaling 
law (0) will be called in the following the A-scaling law. 
We shall return to this form later when discussing the 
second-order critical behavior in finite systems. 

The case < ^ < 1 is more subtle, because < Mx > 
cannot be defined. In this case, a limiting behavior is 
expected to take a form : 



Ni/f^f^ 



f m \ 



(9) 



There is no connection between < Mx > and A^ in (0). 
Instead, one can show that the typical value of Mx for 
large values of N is 

is : 



AI^ ~ N^/f^ , so the correct scaling 



M*xfM^{m)^P^.,p 



M 



N 



(10) 



Notice that |lO[ ) resembles the limit when /z ^ 1 of the 
scaling law (H). This scaling form is not very useful be- 
cause M^ is not as precisely defined as the average value. 
It shows however, that it does not make sense to search 
for the scaling relations such as m) with the scaling ex- 
ponent 1/iJ, larger than 1. 



III. ORDER PARAMETER FLUCTUATIONS 

Fluctuations of the order parameter are expected to 
have different properties at the critical point and out- 
side of it. Far from the critical point, the correlations 
are short-ranged and the results of the previous chapter 
apply in this case. On the contrary, fluctuations of the 
order parameter close to the critical point are correlated 
throughout the whole system. Hence, the 'macroscopic' 
order parameter cannot be written as a simple sum of un- 
correlated block variables, and its asymptotically stable 
distribution differs from those discussed in chapt. 2. In 
the next section we shall discuss some features of these 
distributions and show that they satisfy the A-scaling 
law. 



A. Finite-size scaling argument in the 
thermodynamical systems 

In the thermodynamic limit, the free energy density of 
an equilibrium system close to the critical point scales as 



f{X\Xe)^X^-'^f{rj,e) 



(11) 



where a, (3 are the usual critical exponents, 77 is the in- 
tensive order parameter and A is the scale parameter, 
e is the control parameter which equals at the critical 
point. Even though finite systems do not exhibit the crit- 
ical behavior, nevertheless their properties resemble those 
of infinite systems if the correlation length ^ is larger or 
comparable to the typical length L of the system. In 
this case, one speaks about the pseudocritical point at a 
distance : 

from the true critical point Q. A^ is the size of the D- 
dimensional system and c is some dimensionless con- 
stant which is negative, if a maximum of the finite- 
size susceptibility lies in the ordered phase, or positive, 
if this maximum is in the disordered phase. One can 
then derive the finite-size scaling of the total free energy 
F{ri^ e, TV) = N f{rj, e) at the pseudocritical point : 



In deriving (^2|), we used the hyperscaling relation : 2 — 
a = vD . The canonical probability density of the order 
parameter Pn[v\ is given by : 



PnU = ZN-^eM-PTF{r^,e,N)] 



(13) 



where the coefficient /3t(= ^/T) is independent of 77 (T 
is the temperature of the system). Using eq. (13), one 
may calculate not only the most probable value of the 
order parameter, which is the solution to the equation 
dPN[ri\/dri = 0, but also the average value of the order 
parameter and the partition function 



^N 



N- 



<ri> 



V 



(14) 



■q* denotes the most probable value of the order param- 
eter. < rj > vanishes for large values of N, since both /3 
and 2 — a = 2/3 -|- 7 are positive. The probability density 
Pn Vi\ obeys the first-scaling law : 



with 



< 77 > PtvW = *(??/< ^ >) = ^{z) 



$(z) ~exp(-/3T/(az,c)) 



(15) 



(16) 



The logarithm of scaling function ^(z) corresponds to 
the non-critical free energy density at the renormalized 
distance e = c from the critical point. The scaling law 
([l5|), which is analogous to (|l^), can be rewritten as : 



< 771 > PtvH =*(^(1)) > 



(17) 



where m = Nrj is the extensive order parameter and Z(i) 
is : 



Z(i) = {m — 777.*)/ < m> 



(18) 



m* denotes the most probable value of m. The scal- 
ing domain is defined by the asymptotic behavior of 
Pisi[m] when m -^ 00 and < m >^ 00, but ^(i) has 
a finite value. Assuming that the scaling framework of 
the second-order phase transition holds, the scaling rela- 
tion (p^ ) is valid independently of the explicit reasons of 
changing < 777 > , and independently of any phenomeno- 
logical details. In other words, an explicit relation be- 
tween the size A'' of the system and < m > need not to 
be known. 

If the order parameter is related to the multiplicity of 
fragments or produced particles, like in the fragmenta- 
tion - inactivation - binary (FIB) process |]10[, then the 
first scaling (|ll) becomes identical to the well known 
KNO-scaling lia]. Only in this case, the multiplicity of 
produced particles is the relevant observable. 



F,(77,iV)-/(7?7V— ) 



(12) 



B. Tail of the scaling function 



C. The A-scaling law^ 



What are the properties of the scahng function for 
large values of mj < m >1 The correct way of ap- 
proaching this problem is to study the system subject 
to a small field h which is conjugate to the order param- 
eter and breaks the symmetry of the distribution. This 
consideration yields [|l3[ : 

<i>(z) -- exp(-az?jt-^) EE exp(-az(^^-|) , (19) 

with 5 = {2 — a — P)/(3. The coefficient a in ( |l9| ) depends 
regularly on the temperature. 

One can express this relation in a different way. The 
anomalous dimension for an extensive quantity m — Nrj 
can be defined as : 



g = lim gjsi = lim 



N 



N^oo d\nN 



(In < m >) 



(20) 



One can see from dlj) and the Rushbrooke relation 
a + 2P + J — 2 , that the anomalous dimension is : 



5 = 1-/3/(7 + 2/5) 



(21) 



Since a and (3 are both positive, therefore g is contained 
between 1/2 and 1 for equilibrium systems at the critical 
point of the second-order phase transition. One may note 
also that : 



1 



1 



9 



P 



>2 



(22) 



Whenever the cluster-size can be correctly defined , like 
for example in the case of percolation or Ising models, 
the exponent r of the power-law cluster-size distribution 
Us ~ s~'^ satisfies : 



1 



T-1 



and 



T- 1 

r-2 



(23) 



The allowed values at the critical point are : 2 < t < 3. 
Consequently, one can learn from the cluster-size distri- 
bution whether the equilibrium system is at the critical 
point and whether the considered extensive quantity can 
be identified with the order parameter. For this kind 
of equilibrium phase transitions , the size of the largest 
cluster is the order parameter since : 



< s„ 



>^7V^^ 



N9 



The cluster multiplicity could be the order parameter 
if T < 2, but this can be satisfied only for certain off- 
equilibrium phase transitions, such as for example the 
shattering phase transition in the FIB process flQ] or in 
the dissipative perturbative QCD |ll[|. 



What happens if the observable m is the TV-dependent 
function of the order parameter : 



77 with ai > g 



(24) 



Writing ( |17| ) with m instead of ry and taking into account 
that PNirjldri = Pat [to] dm, one finds : 



< TO >'^ PatH = $(Z(A)) , 

with the scaling variable : 

■2(A) = (to — TO*)/ < TO > , 

where : 

A = g/ai < 1 . 



(25) 
(26) 
(27) 



The normalization of Pat [to] and the definition of < ttt, > 
provide the two constraints which are consistent with : 
< A < 1. This A-scaling law (0 < A < 1) is satisfied, 
as shown in sect. 2.5, by the asymptotically stable dis- 
tributions with 1 < /i < 2 (eq. (pf)). The asymptotically 
stable distributions P^./3 with < /i < 1 satisfy the first 
scaling law, i.e. the limiting case of the A-scaling for 
A = l. 

The scaling function ^(^(a)) has the identical form as 
<i>(z(i')), except for the inversion of the abscissa axis. In 
particular, its tail for large Z(a) has the same form : 



$(z(A)) ^ exp 



(--(A)) 



exp 



-(A) 



(28) 



as given in (|T^) . Notice that the A-scaling of an extensive 
variable : to — N(l — r]) = Nfj, can be used to determine 
directly the anomalous dimension, since in this case : 
A = g. At the phase transition : < Nfj > ~ iV, but the 
finite-size corrections are algebraic. One should mention 
in passing that if ai < g, then < m >^ N < rj > and 
A = l. 

It is important to realize that the A-scaling (0 < A < 
1) and the first-scaling (A = 1) laws are satisfied by the 
distributions of the quantities which are the sum of ei- 
ther independent or correlated random variables. In this 
sense, these scaling laws are the fundamental character- 
istics of the composed statistical system. 



D. Cumulant relations 



Let us now introduce the cumulants Kq from the for- 
mula for the moment expansion of the probability distri- 
bution P[m] : 



In 2_. P['m] exp{mu) 

\m=0 J 



g=0 ^ 



(29) 



In the case of A-scaling law (eq. (p5[)), the cumulants 
satisfy : 

Kq -< TO >«^=< TO >«^ (g > 2) . (30) 

In the case of the first-scahng law (eq. ([l7|)), one finds : 

Kq r^<m>'i . (31) 

Notice also that the scaled variance k,2/ k\ is independent 
of iV for A = 1 and tends to when TV ^ cx) for < 
A < 1. In this sense, the scaling law ( |l7| ) is the large 
fluctuation limit of the probability distribution Pm [to] . 

E. The intermittent behavior 

Intermittency in particle physics is defined usually as 
the scale-invariance of the factorial moments with respect 
to the change of the size of the phase-space cell 5y : 



FgiSy) (X (Sy)-^^ 



(32) 



when Sy —f 0. The exponent (j)q is called the intermit- 
tency index. For the multifractal distribution we have 



J2pm)=<Pr\Sy)>^{Syy'-'^^ , 



(33) 



where Pi{Sy) is the probability for the particle to be in 
the bin i of size 6y and r(g) is given by : 



Since : 



Tiq) = iq-l)DqA 



Fq^iSy)-^^-'-^^ <pi''> , 



(34) 
(35) 



where D is the support dimension, therefore : 

C^q = {q-l){D-DqA) . (36) 

The quantity I?gA appears here as the apparent Renyi 
dimensions. Hence, knowing the intermittency indices 
and the value of A, one can extract the 'true' multifractal 
spectrum. 

F. Scaling laws and the self-similarity 

The self-similarity of the statistical system which obeys 
the first-scaling law (eqs. (jlTl), (fsf)) means that : 



< Nri >xN^ A^ < Nr] >Ar= X" < Nrj > 



N 



The Hurst exponent is equal to the anomalous dimension 
g (eq. (pOJ)). In the case of the A-scaling, H = g = 
aiA (see eq. (|27|)). The exponent ai is the apparent 
anomalous dimension for m defined by (p4) : < m >ajv^ 



A°i < m >N. It is easy to see, that these properties are 
analogous to those we have discussed in sect. 2.4 with 
Yat corresponding to Nr; and ^ Xj corresponding to to. 
Namely : 

01 = 1, g = l/fJ., A = 1/^ 
for 1 < ^ < 2, whereas 

ai = 5 = 1/m, a = 1 
for < ^ < 1. 

G. The off-critical scaling 

Away from the critical region, if m is the order param- 
eter then the finite system exhibits the second-scaling law 
(the fimit A = 1/2 of the A-scaling law (H), ( p6|)) in the 
ordered phase, and the first-scaling law (eqs. (|17|, ([l^)) 
in the disordered phase. In both cases, the tail of the 
scaling function is gaussian {i> = 2) and the finite-size 
corrections are exponential. The second-scaling law has 
been found in the shattering phase of the non-equilibrium 
FIB process flj] and in the 'liquid' phase of the equilib- 
rium percolation process p5| . The first-scaling law with 
the gaussian tail has been found in the disordered phase 
of the Ising model mm . 

One would like to know how a finite system evolves 
close to the critical point when the control parameter e 
tends slowly to : e - 7V2ai-2 ^^ ^ -^-^ ^ j^^ ^-j^jg ^^^^ 

< m > and m* have a common behavior ~ N"''^ , and the 
probability distribution Pat [to,] satisfies the scaling : 



< m >^ PnViA ~ exp 



(to. — m*y 



< m > 



2A 



(37) 



with c > and A = 3ai/2 — 1 . The case ai — 1, i.e., e = 
const, corresponds to the second-scaling law. ai = 3/4 
corresponds to the first-scaling law since the finite system 
is in the critical region (oi ~ g). In between these two 
limiting cases, the A-scaling holds with 1/2 < A < 1. In 
all cases, the scaling function has a gaussian form. 



H. Equilibrium vs off-equilibrium critical systems 

Let us consider the bond percolation model. In the 
regular lattice, each site corresponds to a monomer and 
a proportion p of active bonds between occupied sites 
is set randomly. Such a network results in a distribu- 
tion of clusters which are defined as the ensemble of 
occupied sites connected by active bonds. For a criti- 
cal value of active bonds Pcr, an 'infinite' cluster (the 
gel) spans the whole lattice. Infinite in this context 
means that the gel contains a finite fraction of the to- 
tal mass of the system. The sol - gel phase transition can 



be suitably studied using moments of the number-size- 
distribution Us : Mi = yi„^„ s^Us , where the sum 
runs over all clusters with the exception of the largest 
cluster s = Smax and the superscript ' recalls this con- 
straint. The mass of the largest cluster is then : N — M[. 
In the infinite system, one usually works with the normal- 
ized moments of the concentration-size-distribution c, : 



E. 



s'^C, 



tj. — ^g^g iT i-s • The concentrations are normalized 
such that : Cg = liniAr^oo rig/N . The probability that a 
monomer belongs to the infinite cluster (the gel) is equal 
to l — m'i with : mj. = limjv^oo M^/N . Therefore, in the 



thermodynamic limit 



1 for p < Per and m'l < 1 



for p > pc- , i.e., a finite fraction of the total number of 
sites belongs to the infinite cluster. 

After this introduction, let us discuss the bond- 
percolation on the Bethe lattice with a coordination num- 
ber z |17| . The concentration-size-distribution in this 
model is : 

cs-z ^^!^";,^^"l,, p'-'il-p)^'-'^'^' (38) 

{{z — 2)s + 2)!s! 

and the first normalized moment equals : 

I-P 
1-p* 

with p* being the smallest solution of the equation : 

p*{l-p*f-'^p{l-pf-' . (39) 

For p < Per — (z — 1)^^, the only solution of eq. ( p9| ) is : 
p* — p, but when p is larger than pcr, then there exists a 
smaller non-trivial solution which behaves as Pcr — \p—Pcr\ 
near Per- Above this threshold, the moment m'l is smaller 



than 1 and behaves approximately as 



l-2(p- 



Per) / {'i^ - Per) ■ 

For large values of the size s, the concentrations 
can be rewritten using the Stirling approximation : 



-5/2 



exp(— as) 



(40) 



a in (EO) equals : 



a = In 


Per \1-Perj 


~(z 


- 2) ln(l - p) < . 



(41) 



Hence 



-5/2 



a-p) 



-s{z-2) 



At the percolation threshold, the concentrations follow 
the power-law and outside this threshold, an exponential 
cut-off is always present. 

If we impose that Smax < oo, even in the thermo- 
dynamical limit, then the activation p corresponding to 
Smax is smaller than pcr and the limiting activation po 



equals Per- The limiting activation pq in percolation is 
analogous to the limiting temperature Tq in the statisti- 
cal bootstrap model (SBM) ^. Both po and Tq follow 
from the consistency conditions in the two models and 
describe the asymptotic mass spectra. Obviously, the 
limiting value Po/To of the control parameter p/T does 
not have to correspond to the critical point in the equi- 
librium system. For this to be true, the exponent r of 
the mass spectrum must be contained in between 2 and 
3. 

One needs two critical exponents to describe com- 
pletely the critical features. In percolation these are for 
example r and a (the exponent of the mean cluster-size 
divergence). The A-scaling analysis allows to determine 
only this exponent which is related to the anomalous di- 
mension (e.g. T in the percolation problem). 

The insights gained from the numerical simulations of 
non-reversible aggregation equations (the Smoluchowski 
model) and FIB equations, both describing the off- 
equilibrium phenomena, and the evidences from analyti- 
cal solutions for gelling and non-gelling systems iQ , pro- 
vide strong indication that the scaling features, which 
were discussed in chapt. 2 for the distribution of the sum 
of random short-ranged correlated variables and in chapt. 
3 for the distribution of the long-ranged correlated vari- 
ables at critical point in equilibrium systems, are valid 
also for non-equilibrium systems. This universal behav- 
ior of equilibrium and non-equilibrium systems, has the 
deep foundation in the relation between renormalization 
group ideas for self-similar systems and the limit the- 
orems of the probability theory for the A-scaling laws, 
both for correlated and uncorrelated variables. The con- 
cept of statistical equilibrium does not intervene at this 
level. 



IV. WHERE WE ARE, WHERE DO WE GO? 

There are two generic classes of dynamical critical phe- 
nomena, which are characterized by different 'relevant 
observables'. The first one is the sequential cluster frag- 
mentation, where the average cluster size decreases and 
the cluster multiplicity increases during the process. The 
cluster (particle) multiplicity is an order parameter for 
this class of critical phenomena. The second one is the 
sequential cluster aggregation, where the average cluster 
size increases and the cluster multiplicity decreases dur- 
ing the process. The size of the largest cluster is here 
an order parameter. In the first class, one finds for ex- 
ample the shattering phase transition in the FIB pro- 
cess flQ] as well as in the dissipative pcrturbative QCD 
process |ll|. The second class contains for example the 
liquid-gas phase transition (Fisher model), the percola- 
tion transition or the sol-gel phase transition in the irre- 
versible kinetic aggregation. Which of them is approxi- 
mately realized in the multiparticle production process, 



can be discovered by studying the A-scaling law, the tail 
of the scaling function, and the anomalous dimension, as 
described in the previous chapters. 



A. The perturbative QCD equations with the 
inactivation mechanism 

Let us call Pn [n; t] the probability to get a cluster mul- 
tiplicity n at time t, starting from an initial cluster of size 
A^ at i = 0. The time evolution equation for the multi- 
plicity is given by the non-linear rate equation llOl : 



dG 



N 



dt 



N-l 

u,t) ^Yl ^J,N~j[Gj{u,t)GN-jiu,t) - 
- GN{u,t)] +In[1 + u~ GNiu,t)] 



(42) 



with the initial condition : Gi{u,t) = 1 + u, and the 
normalization condition : Gn{u — 0,t) — 1. Gj\j{u,t) is 
the generating function of the probability distribution : 
GN(u,t) = J2^=i-PN[n;t]{l + uy^. The partial derivative 
in (^2|) is taken at a fixed size N. The sum on the right 
hand side (rhs) of eq. (|4^) represents a binary fragmen- 
tation of the primary cluster N into the daughter clusters 
of mass j and N — j respectively. The second term on the 
rhs side is the dissipative term which is responsible for the 
inactivation. One can transform the discrete variable j in 
( [42| ) into a continuous one z = j /N, which varies from 
to 1. With this change, the fragmentation kernel J^j^n-j 
is replaced by a splitting function : '^z.i-z = ^^j-N-j- 

In the limit of 'no dissipation', the FIB process yields 
the rate equations of gluodynamics in the Next to Next to 
Leading Logarithm Approximation (NNLLA) jlj]. The 
inactivation term of FIB model yields a unique prescrip- 
tion of how to obtain the rate equations of perturbative 
QCD in the NNLLA, which would contain the hadroniza- 
tion effects phenomenologically through the inactivation 
mechanism of parton cascading. 

The time t in eq. (p3) arises within the fragmentation 
and inactivation kernels , which themselves are probabili- 
ties per unit of i. We define then the time as : i = TlnF, 
where T is a constant, Y — ln(A^6/Qo) , Qo =const and 
Q plays the role of time and orders the sequence of events. 
Assuming now that all physical quantities depend only 
on the variable Y and not on N and O separately, we 
transform (E2f) into : 

X GiY + log(l - z),u) - G{Y, u)]dz 

+ niY,u) , (43) 



where : 70^(1") = 2Tr-'^Ncas{Y) ^ T/Y, and OsiY) is 
the QCD running coupling constant. I{Y) in (Q, is 
the inactivation function written in the new variables. It 
should be stressed that no cutoff is needed and the par- 
ton cascades are inactivated gradually by the inactivation 
function T{Y). The generalization of eqs. (Esl) to include 
both quarks and gluons can be found in |l] 



B. Fragmentation scenario for e^e collisions? 



with 



7^(r, u) = joHy)Ay) [1 + u- G(r, u)] 



(44) 



Eqs. {h^ of the dissipative gluodynamics have been 
applied for the description of the hadron multiplicity data 
in e+e" collisions [|ll|. The available data have been well 
described assuming the gaussian inactivation mechanism 
of the parton cascades. Based on that, it was concluded 
that the multihadron production in e~^e~ collisions re- 
sembles the critical off-equilibrium shattering process and 
the non-perturbativc hadronization stage, approximated 
by the inactivation function I{Y), is slightly shifting this 
process from the critical line into the shattering phase. 
The small deviations from the first-scaling (the KNO- 
scaling in this case) seen in the data for y/s < 100 GeV 
is a finite-size effect related to the small size of an ini- 
tial jet. The asymptotic region could not be precisely 
determined in these studies but, nevertheless, the limits 
imposed by the experimental data on the parameters of 
the inactivation function provide the upper limit for the 
borderhne of an asymptotic region : y^ ^10^ TeV. 

The experimental analysis of the ratio of the aver- 
age charged multiplicity to the dispersion in e~^e~ —^ qq 
events shows that the multiplicity distributions obey the 
first-scaling law [pO| . This confirms the relevance of the 
hadron multiplicity observable in this process. Unfortu- 
nately, the tail of the multiplicity probability distribu- 
tion has not been resolved experimentally, so it is un- 
known at present whether e~^e~ — *■ hadrons is the crit- 
ical process at presently available energies and whether 
the hadron multiplicity is the order parameter in this 
process. Nevertheless, even though the experimental 
data are yet incomplete in many respects, the fact that 
such a simple observable as the hadron multiplicity is 
relevant allows to hope that the future studies at still 
higher collision energies will solve the enigma of pos- 
sible off-equilibrium shattering phase transition [lC| ] in 
the process : e"'"e~ — > hadrons. This specific transi- 
tion, dominated by the perturbative regime of the QCD, 
is extremely interesting even though unrelated with the 
quark-gluon phase transition, which in turn is dominated 
by the strong-couphng regime of the QCD. 



C. Aggregation scenario for pp and A A collisions? 

Hadron-hadron and nucleus-nucleus collisions show 
quite distinct features from those of lepton-lepton col- 



lisions. First of all, the multiplicity distributions show 
a clear deviation from the first-scaling law [£1|. The 
pp data for multiplicity distributions at y/s — 200, 546 
and 900 GeV Q can be approximately superposed at 
around the maximum in the variables of the A-scaling 
law with A = 0.9 122]. On the other hand, the strong 
scaling violations are seen in the tail for large multiplic- 
ities showing that the global scaling of multiplicity dis- 
tributions for different ^fa may not be possible. In this 
context, the two-component scenario including soft and 
semi- hard classes of events has been discussed Q . The 
system size-dependence of these components cannot be 
accurately determined based on the existing data, so it 
is at present difficult to establish a physical meaning of 
this scenario. One should however mention, that finite 
systems in the neighborhood of the first-order liquid-gas 
phase transition, show a similar two-component behavior 
in the cluster multiplicity. 

Correlation studies have disclosed the deficiency of the 
Lund Monte Carlo program PYTHIA containing jets, 
resonance production and the Bose-Einstein correlations 
as final state interactions, for the description of the differ- 
ential density correlation function p4| . The correlations 
are both overestimated at high -pT (high multiplicity) and 
underestimated at low p^ (low multiplicity) by PYTHIA 
p4| , p5[ . Significant improvement can be gained assum- 
ing existence of clusters of low pT |2^,^,^ , as included 
in the non-diffractive event generator GENCL [E6[. The 
picture of hadron-hadron collisions which emerges from 
these phenomenological studies shows the necessity of 
a combined description of hard scattering processes fol- 
lowed by the string fragmentation (the perturbative QCD 
regime) , and the formation and decay of clusters in the 
non-perturbative QCD regime. The absence of an ade- 
quate computer program simulating these two aspects of 
hadron-hadron (nucleus-nucleus) collisions is presently a 
principle obstacle, slowing down a progress in the under- 
standing of the production mechanism of soft hadrons in 
the strong-coupling, long-distance regime of QCD. The 
theoretical activity aiming at the development of such 
program should be accompanied by experimental studies 
of the appropriate methods to evidence the existence of 
clusters and to study their sizes. 

The arguments put forward in favour of the cluster 
picture refer often to the SBM jl^. In this model, the 
number of species of possible constituents of the fireballs 
and the number of species of fireballs grow asymptoti- 
cally like ^ 771"'^ exp (to/To )i with r = 3 [ESl. As dis- 
cussed in chapt. 3, this value of r excludes a possible 
identification of Tq with the critical temperature of the 
thermodynamical quark-giuon phase transition. Another 
family of models refers to the spontaneous breakdown of 
the chiral symmetry of the QCD vacuum at a critical tem- 
perature. In the approximation to„ k, m^ w 0, rua — > cx), 
the chiral QCD phase transition is of second order and 
the order parameter is given in terms of the sigma and 



pion field condensates |2^. In general, we get the guid- 
ance in designing appropriate models simulating the crit- 
ical behavior in finite temperature QCD from the lattice 
QCD calculations (see e.g. fS^]). In all these models, the 
proper characterization of the soft phase of QCD should 
involve an understanding of features of the largest cluster 
(fireball) distribution. 

To design the relevant experimental observable in pp 
and AA collisions is by far the most difficult task. At 
the same time, this objective is the most important chal- 
lenge of soft physics which if unsolved will jeopardies the 
analysis of different phases of hot hadronic and quark- 
gluon matter. The hadron (baryon) multiplicity distri- 
bution may contain useful information only if the studied 
system happens to be found at the critical point, what 
obviously limits the relevance of this observable. In this 
case, the the study of the scale dependence of correlations 
by means of the scale dependence of factorial (cumulant) 
moments of the density distribution, could reveal the un- 
derlying fractal structure of the largest cluster |Q . This 
is for the moment largely a theoretical possibility be- 
cause the exact reconstruction of the geometry of the 
largest cluster may be even a harder experimental task 
than finding its mass event by event P4] . 



V. CONCLUSIONS 

There is still a long way to understand the 
multihadron-production processes and, in particular, 
those aspects of these processes which are related to 
the strong-coupling long-distance regime of QCD. Ex- 
perimental studies of correlations in small domains of 
the phase-space, posed the problem of the information 
contained in the multiparticle fluctuations of produced 
hadrons and, more generally, in the fluctuations of ob- 
servables. We have shown that these fluctuations have 
universal features even in small finite systems, indepen- 
dently of whether the studied process is classical or quan- 
tum, equilibrium or non-equilibrium, dissipative etc. The 
information contained in the particle density and corre- 
lations can be extracted by determining (i) the scaling 
features (the A-scaling law) of the measured probability 
distributions of the observable, in particular, the asymp- 
totic properties of the scaling function, (ii) the anomalous 
dimension, and (iii) the scaling features of multiparticle 
correlations (clustering) in small phase-space cells. In 
this way, the relevance of the observable can be tested 
and/or the constraints on the multihadron-production 
mechanism can be determined. 

Much of the experimental search for the scaling fea- 
tures in the multiparticle distributions has been moti- 
vated by the prediction of the KNO-scaling as an ulti- 
mate symmetry of the S'-matrix in ultrarelativistic colli- 
sions. This scaling is however only a particular case of 
the general A-scaling law which has been introduced only 



recently. For that reason, neither a systematic analysis 
of the deviation from the KNO-scaling in terms of the 
A-scaling, nor the asymptotics of the scaling functions 
are available at present. 

A crucial problem is the determination of the rel- 
evant observable for each multihadron-production pro- 
cess, which could exhibit different phases of the studied 
process through the non-trivial fluctuation pattern and 
its evolution. Whereas the hadron multiplicity distri- 
bution seems to be relevant for e~^e~ collisions, which 
share many aspects of the fragmentation scenario, this 
observable seems to be of secondary importance in pp 
and AA collisions. In these latter collision processes, 
it is plausible that the aggregation scenario dominates 
and, hence, the cluster measures of the data have to be 
studied. The development of new methods of extraction 
of the 'biggest cluster' event by event, accompanied by 
the development of reliable programs simulating both the 
perturbative QCD regime and the clustering aspects in 
the non-perturbative regime, are the challenging, urgent 
problems which will determine the future evolution of the 
soft physics and the multiparticle correlations. 
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